Absence of Wigner molecules in one-dimensional few-fermion systems with 

short-range interactions 



Jing-Jing Wang/ Wei Li,^ Shu Chen,^ Gao Xianlong,^'!^] Massimo Rontani,^'[^ and Marco Polini^ 

^Department of Physics, Zhejiang Normal University, Jinhua, Zhejiang Province, 321004, China 
^Department of Physics, Fudan University, Shanghai, 200433, China 
^Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
^S3, Istituto Nanoscienze-CNR, Via Campi 213/A, 41125 Modena, Italy 
^NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, 56126 Pisa, Italy 

We study by means of exact-diagonalization techniques the ground state of a few-fermion system 
with strong short-range repulsive interactions trapped by a harmonic potential in one spatial dimen- 
sion. Even when the ground-state density profile displays at strong coupling very well pronounced 
Friedel oscillations with a '4/cf periodicity', the pair correlation function does not show any signa- 
ture of Wigner-molecule-type correlations. For the sake of comparison, we present also numerical 
results for few-electron systems with Coulomb interactions, demonstrating that their ground state 
at strong coupling is, on the contrary, a Wigner molecule. 
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Introduction. — Understanding the nature of the 
ground state of strongly correlated few-fermion systems 
has attracted a great deal of interest. A large body 
of literature is available in the context of few-electron 
quantum dots, or 'artificial atoms' pLj. In this case the 
relevant interaction potential is the long-range Coulomb 
force between electrons. It is by now well established 
that when Coulomb interactions are sufficiently strong 
the ground-state of a quantum dot — either one- [2 or 
two-dimensional [3 — is a 'Wigner molecule' (WM). This 
jargon stems from the Wigner crystalline order that is 
displayed by an electron gas at ultra-low densities [4]. 
The WM state is not characterized by a spontaneous 
breaking of translational symmetry (which is absent to 
begin with due to the confining potential of the dot) but 
rather by strong short-range order that is evident in the 
density-density correlation function (the so-called pair 
correlation function). The existence of WMs has been 
recently confirmed by inelastic light scattering studies of 
few-electron circular quantum dots [5 as well as by tun- 
neling spectroscopies of quantum wires [6] and carbon 
nanotubes [7]. 

The ground state of one-dimensional (ID) systems of 
fermions trapped in a parabolic potential and interacting 
through a short-range potential has also been the subject 
of numerous analytical and numerical studies [H 1^ . The 
interest in these systems is not merely academic since 
gases of fermionic atoms trapped in 'atomic quantum 
wires' can be realized experimentally [10 and cooled to 
temperatures T ~ 0.1 Tp, where Tp is the Fermi tem- 
perature [11 . The relevant interaction potential in this 
case is a zero-range Fermi pseudopotential acting only 
between antiparallel-spin fermions [12], V{x) = gm5{x)^ 
with an amplitude which is controlled by the 
three-dimensional (3D) scattering length. Recent ad- 
vances in atomic physics have made possible also to con- 



fine only a few fermionic atoms with tunable short-range 
interactions in optical dipole traps [14 . 

In the non-interacting limit the ground-state 

density profile of TV fermions trapped in a parabolic po- 
tential displays, despite the smooth boundaries, 'Friedel 
oscillations' with N/2 peaks [15], which are ultimately 
due to the Pauli principle. In the non-interacting 
limit, indeed, the ground-state density no(x) is ob- 
tained by occupying exactly N/2 harmonic-oscillator 
levels, each one with two fermions with antiparallel 
spin: no(x) = 2 ^^^^""^ |(/)n(x)p, where (/)n(x) = 
(2^n!7ri/2^ho)"^/^ exp (-^V2) i^n(0 with ^ = x/4o and 
4o = {h/mujY^'^ are the eigenfunctions of a ID harmonic 
oscillator with frequency uo. Here Hn{C) is a Hermite 
polynomial of degree n. Straightforward mathematical 
manipulations yield the following approximate expression 
for no(x) away from the trap edges [16^: 
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where nTF{x) = (2Ltf/C)[1 - (^/^tf)^]^/^ is the 
Thomas- Fermi density profile, Ltf = 4o^^^^5 and 
^f(O) = 7rnTF(0)/2 = N^^'^/iho (the expression for kY{x) 
away from the trap center can be found in Ref . [T6]) . 

What happens to this simple single-particle physics 
when the strength of inter-particle interactions is in- 
creased? In the limit giB oo the delta-function in- 
teraction imposes an effective Pauli principle between 
antiparallel-spin fermions. In this limit we thus ex- 
pect oscillations in the ground-state density noo{x) with 
a bulk periodicity controlled by 4/cf(0). These have 
been dubbed in the literature 'Wigner oscillations'. The 
transition from the 2A:f(0) Friedel oscillations to the 
4/cf(0) Wigner oscillations in a parabolic trap is a smooth 
crossover: see, for example, the extensive density-matrix 
renormalization group study by Soffing et al. 
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Looking, however, at a one-body observable such as the 
ground-state density profile does not shed light on the na- 
ture of the ground state at strong coupling. Especially 
important in low-dimensional systems are indeed correla- 
tion functions. In the absence of the parabolic trapping, 
we know from bosonization [17] that at strong coupling, 
i.e. when Kp < 1/3, Kp being the Luttinger liquid pa- 
rameter in the charge sector, ID interacting fermions are 
dominated by 4/cf charge-density- wave (CDW) correla- 
tions. Since in a Galilean invariant system Kp = v-p/vp^ 
where is the bare Fermi velocity and Vp is the ve- 
locity of density excitations, and Vp 2vy in the limit 
giB oc [8j, we expect that the ground state at strong 
coupling is completely dominated by 2/cf spin-density- 
wave correlations with no space for incipient CDWs. The 
situation changes dramatically in a system with long- 
range interactions, which stabilize 4/cf CDW correla- 
tions [IMS]. 

In this Letter we analyze the nature of the ground state 
of a ID few-fermion system at strong coupling in the 
presence of a parabolic potential, which breaks Galilean 
invariance. More precisely, inspired by the body of litera- 
ture discussed above, we study both ground-state density 
profiles and pair correlation functions. We focus on the 
strong coupling regime, which is difficult to access with 
analytical techniques, and rely on accurate numerical cal- 
culations based on the exact-diagonalization method (aka 
'full configuration interaction') [2(1. We find that even 
when the ground-state density profile displays well pro- 
nounced Wigner oscillations with a 4/cf(0) periodicity, 
the pair correlation function does not show any signa- 
ture of WM-type correlations. We then highlight the 
role of the interaction range by comparing these findings 
with those for a system of electrons interacting through 
the Coulomb potential. We demonstrate indeed that the 
ground state of this system at strong coupling is a WM. 

The model and the observables of interest. — We con- 
sider a two-component Fermi gas with N atoms con- 
fined inside a strongly elongated harmonic trap. The 
two species of fermionic atoms are assumed to have 
the same mass m and different pseudospin a =t or | 
( hyp erfine- state label). The trapping potential is axi- 
ally symmetric and characterized by angular frequencies 
00 ± and uj in the radial and longitudinal directions, re- 
spectively, with UJ <^uj±. Correspondingly, we introduce 
the harmonic-oscillator lengths a_\_ = {h/muj±_Y^'^ and 

4o = {n/mujf/\ 

The gas is dynamically ID if the anisotropy param- 
eter of the trap is much smaller than the inverse atom 
number, cj/cj^ <C N~'^. It can thus be described by an 
inhomogeneous Gaudin-Yang Hamiltonian, 

fc2 r+oo 

?i = t + V + W = -^^y_ dx¥,{x)dli^,{x) 
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where ^J-(x) [^cr(^)] is a field operator that creates (de- 
stroys) a fermion with spin a at position x and ^id — 
Ah^a^c/ {ma\) (in the limit age <C a^) is a parameter that 
determines the strength of inter-particle repulsions [13]. 
The 3D scattering length age can be tuned by means of a 
magnetic field [10] . The first term in Eq. ^ (T) is the ki- 
netic energy whereas the second term (V) describes two- 
body short-range interactions. Finally, the third term 
(W) is the parabolic trapping potential. 

Choosing as unit of length and hw as unit of en- 
ergy, the Hamiltonian Q is governed by the dimension- 
less coupling parameter 



A: 
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In this Letter we focus our attention on the dependence 
of two key quantities on A for strong inter-particle re- 
pulsions (A ^ 1): i) the local spin-resolved density, 
ncT{x) = (^J.(x)^cr(^))/^cr, and ii) the pair correlation 
function (PCF) 



+ 0O 



dx' Qaa'ix' +2'^' 
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which is defined in terms of the two-body correlator 



:{x,x') 
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Ncr being the number of atoms with spin a. We also intro- 
duce the total density, n{x) = Ncrncr{x) /N ^ and the 
total PCF, -g{x) = E.,.^ N,{N,> - S^^^)g^^^{x)/N{N - 
1). In all the definitions above (...) denotes the expec- 
tation value over the ground state of the Hamiltonian 
([2|. We remind that gcra'{x^x^) measures the conditional 
probability of finding a fermion with spin at position 
x^ when another fermion of spin a is known to be at 
position X, whereas gaa'{x) provides the probability of 
finding two fermions at a relative distance x. 

For the sake of comparison, we present also numerical 
results for a few-electron system in the same trapping 
potential. In this case the second quantized Hamiltonian 
is identical to the one in Eq. ^ except for the term in 
the second line, which now reads: 



1 p-\-oo r-\-oo 

= oE / / dx' ¥,ix)^Ux') 

a, a 

X Vc{\x-x'\)^^,{x')^,{x) . (6) 

As a model long-range potential we use the following reg- 
ularized Coulomb interaction 



Vc{\x-x'\) = 



K^a^ + {x — x'Y 



(7) 
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where a is a real positive quantity with physical dimen- 
sions of length and k is the relative dielectric constant 
of the host semiconductor. The Hamiltonian He = 
T + Vc + W describes few-electron quantum dots em- 
bedded e.g. in a ID quantum wire [2 . Using the same 
dimensionless variables which have been used above for 
the inhomogeneous Gaudin-Yang Hamiltonian, we find 
that the dimensionless coupling constant that determines 
the strength of electron-electron interactions is given by 



Ac 



(8) 



Due to the presence of the parameter a in Eq. ([7|), the 
few-electron problem does not depend only on Ac but 
also on the ratio a/£ho, whose value has been fixed for all 
our numerical calculations to a/^ho = 10~^ [21j. 

Numerical results. — In Fig. [T] we report numerical re- 
sults for the spin-resolved density profiles n(j{x) oi N = A 
and = 5 fermionic atoms. All the data presented in 
this figure and in all the figures below have been calcu- 
lated for A = 15 and are qualitatively similar to those 
obtained for other values of N ^ whereas A/" > 6 is be- 
yond our reach. In this ultra-strong coupling regime 
mean-field-like methods yield completely incorrect re- 
sults (while we have checked that our numerical results 
at weak coupling can be explained very well by Bethe- 
Ansatz density-functional theory P). In both panels of 
Fig. [l] we clearly see N distinct peaks in the total den- 
sity profile n{x)^ i.e., Wigner oscillations with a 4/cf(0) 
periodicity. Note also the 'antiferromagnetic' spatial pat- 
tern of the spin-resolved densities n^^x) in which n^{x) 
has a maximum in correspondence of every minimum of 
n^{x). Quite surprisingly, the PCFs corresponding to 
the density profiles in Fig. [l] which are plotted in Fig. |2j 
do not display any sign of WM-type short-range order. 
Despite the well-defined Wigner oscillations in the den- 
sity profile, the ground state seems to display liquid- type 
correlations, without any sign of interaction-induced lo- 
calization. 

In Fig. [3] we compare the results reported in Fig.[l]with 
those for few electrons interacting through the Coulomb 
potential. For the sake of comparison, in producing the 
data for Fig. |3] we have chosen Ac = 15, which is the 
same value used for the short-range coupling A in Fig. [T] 
Also the total density profile in Fig. [3] is characterized by 
Wigner oscillations with a contrast though which is much 
higher than in the case of fermions interacting through 
short-range forces. Moreover, at this large value of Ac, 
the ground state of the few-electron system is definitely 
a WM. This is highlighted in Fig. [4] where we have pre- 
sented the total PCFs corresponding to the density pro- 
files in Fig. [3] In the Coulomb-coupled system we note an 
'excluded volume' region, i.e., a finite range of values of x 
close to the origin in which g{x) ~ 0, which is typical of 
electronic systems at ultrastrong coupling [22 . Note also 
that g{x) displays exactly N — 1 high-contrast peaks, a 





FIG. 1: (Color online) Spin-resolved and total ground-state 
density profiles (in units of ^^0^) as functions of x/^ho for A = 
15. Panel a) Results for = 4 particles (all curves overlap). 
Panel b) Results for AT = 5. 



clearcut signature of the WM nature of the ground state. 

In conclusion, we have studied density profiles and pair 
correlation functions of few-fermion systems with short- 
range interactions at strong coupling. We have discov- 
ered that the nature of the ground state of these systems, 
at least in the regime of coupling constants we have inves- 
tigated, is not a Wigner molecule, despite the well defined 
'Wigner oscillations' displayed by the density profile. We 
believe that the ground state of this system is not 'local- 
ized' into a Wigner molecule because of the short-range 
nature of the inter-atom interactions. Systems with iden- 
tical statistics but long-range Coulomb forces, such as 
electrons trapped in a quantum dot embedded in a thin 
quantum wire, are Wigner molecules at similar coupling 
constants. It will be useful to study in the future the in- 
termediate case of dipolar interactions, which are attract- 



4 



1 1 


up up 




up down 




-N. total 




\ \\ Ar = 4. 




a) 


/ 1 1 


1 1 



Short — range interactions 
Coulomb interactions 



3 

x/4o 





0.30 r- 
0.25 - 
0.20 - 
0.15 - 
0.10 - 
0.05 - 
0.00 — 



b) 



Short — range interactions 
Coulomb interactions 




FIG. 2: (Color online) Pair correlation function Qaa'ix) (in 
units of i^^) vs x/iho for the same system parameters as in 
Fig. [l] Note that ^tt(^ = 0) = 0, a manifestation of Pauli's 
exclusion principle. 

ing so much interest in the cold-atom hterature [23ti25] . 
Since one-dimensional few-fermion systems with tunable 
short-range interactions have been recently created [14], 
our predictions can be tested experimentally by using lo- 
cal probes that access both spatial density distribution 
and density-density correlation functions. 
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FIG. 3: (Color online) The total density profile n{x) of a 
system of fermions interacting with short-range interactions is 
compared with that of electrons interacting through Coulomb 
interactions. The data for short-range interactions have been 
produced by setting A = 15 (same data shown in Fig.[T]). The 
data for Coulomb interactions have been produced by setting 
Ac — 15 and a/iho — 10~^. 
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